Expansion of a dipolar condensate 
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We discuss the expansion of an initially trapped dipolar condensate. When the nominal isotropic s- 
wave interaction strength becomes tunable through a Feshbach resonance, anisotropic dipolar effects 
are shown to be detectable under current experimental conditions [E. A. Donley et al., Nature 412, 
295 (2001)]. 
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Atomic Bose-Einstein condensates (BEC) are realized 
at extremely low temperatures jj], @, || , when short range 
atom-atom interaction can be described by a single pa- 
rameter: the s-wave scattering length a sc . This is a valid 
approximation as all higher order partial wave collisions 
die away in a short range potential when the collision en- 
ergy approaches zero. In more realistic models, however, 
even ground state atoms may be polarized such as in an 
external magnetic trap, when the direction of its spin (of 
the valence electron for alkali) becomes aligned. BEC 
with dipole-dipole interaction has attracted considerable 
attention in recent years ||, [|, ^], [?], |8|, |L and trapped 
fermionic dipoles have also been studied [10|. 

Despite the fact that dipolar interactions are now 
widely known to exist in a trapped BEC, it has not been 
directly detected yet. This is because the dipole inter- 
action is much weaker than the contact interaction un- 
der most circumstances. Recently we studied the small 
amplitude shape oscillation of a dipolar condensate 0. 
When the s-wave scattering length a sc is tuned near to 
zero as realizable in the 85 Rb experiment |ll], [l^], we 
showed that the anisotropic dipolar effect becomes de- 
tectable under current experimental conditions. 

The free expansion of an interacting condensate after 
its sudden release from the trap is now a standard di- 
agnostic tool in BEC physics |l3| 111]. It has also been 
extended to the case of an interacting Fermi gas Jl5| . In 
this paper, we investigate the free expansion of a dipo- 
lar condensate jl0|. Our results show that this may also 
lead to the experimental detection of dipolar interactions. 
This paper is organized as follows; First we briefly re- 
view the formulation for a dipolar BEC; we then discuss 
our study of its expansion with both variational and nu- 
merical methods. We will subsequently concentrate on a 
detailed study for situations corresponding to the 85 Rb 
experiment [ |TT| . Then we conclude with some new results 
on the shape and stability of a dipolar condensate. 

For simplicity, we study a trapped dipolar BEC as- 
suming all atomic dipole moments are equal and aligned 
along the z-axis. The atom-atom interaction is then sim- 
plified to be 



V(R) = g S(R) + . 92 
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where R = f — r Or is the polar angle of R, go — 
4irh 2 a sc /M representing the contact interaction, and in 



atomic units <?2 = a 2 (0)£ 2 (a(0) atomic polarizability) 
or [i 2 (/i magnetic dipole moment). The corresponding 
Gross-Pitaevskii equation (in adimensional form) is then 
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where tp(r) is the condensate wave function (normal- 
ized to 1) and a cylindrical symmetric harmonic trap in 
dimensionless unit V^ xt (r) = (x 2 + y 2 + A 2 z 2 )/2, with 
an aspect ratio A. The length unit is a± 
and energies are measured in units of huj± is the 
radial frequency of the trap). Pa = V27rNa sc /a± 
measures the contact interaction strength while P2 = 
\/27rNg2/(3hujj_a 3 L ) denotes the strength of the dipole- 
dipole interaction. The ground state wave function of a 
dipolar condensate can be found by replacing the lhs of 
Eq. (Q) with —flip, where [i is the chemical potential. 

To proceed with the study of the condensate free ex- 
pansion; we can solve the Eq. (||) numerically. We ini- 
tialize the wave function to the self-consistently solved 
ground state in the presence of the trap V ex t, then find 
the time evolved wavefunction from the Eq. (^) by em- 
ploying a self-adaptive Runge-Kutta method without the 
V eK t- In practice, this becomes an expansive calculation 
as we we have to use a rather large spatial grid to accom- 
modate the ever-expanding wave function and to obtain 
sufficient accuracy. We find that it is necessary to check 
the accuracy of solutions repeatedly on different sized 
grids. 

Alternatively, the variational approach as developed 
earlier |t]] can be used. In this case, we assume that the 
wave function always takes the form of a Gaussian, and 
transform the Equation (^J) into the following equations 
for variational parameters 



q r + q r = 
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which in fact represent the radial and axial condensate 
widths of the expanding condensate as q r = W (x 2 ) /2 = 
■J (y 2 )/2 and q z = (z 2 ) /2. k = q r /q z is the condensate 
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FIG. 1: Comparisons of free expansion results from the nu- 
merical (solid line) and the variational (dashed line) ap- 
proaches as discussed in the text, for A = 2, Po = 10, and 
P2 = 10. The upper panel is for the condensate widths (in 
units of a±) and the lower panel for the various energy com- 
ponents (in units of hu±). 



aspect ratio, generally different from the trap aspect ra- 
tio A, and /(«) = [-4k 4 -7k 2 + 2 + 9k 4 P(k)]/2(k 2 -1) 2 , 
g( K ) = [-2k 4 + 10k 2 + 1 - 9k 2 H(k)]/(k 2 - l) 2 with 
H(k) = tanh _1 \/l — k 2 /V1 — k 2 - The kinetic and inter- 
action energy per atom can then be expressed in terms 
of condensate widths according to 
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The validity of the variational solution has been 
checked thoroughly for the frequencies of condensate 
small amplitude shape oscillations || We have 

checked over a wide range of parameters before conclud- 
ing that it is also justified for use in the free expansion 
problem. This enabled us tremendous freedom in ex- 
ploring the expansion dynamics without resorting to the 
time consuming numerical solutions. Figure^ shows the 
time dependence of condensate widths and energy com- 
ponents assuming A = 2 and Po = P2 = 10. We see 
that the variational approach indeed gives a very good 
approximation to the numerical calculations. This result 
is also confirmed for A from 0.1 to 3. 

The time dependent behavior of the energy compo- 
nents is rather interesting. As shown in Fig. 0, irrespec- 
tive of the initial signs of the overall dipole interaction 
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FIG. 2: The time dependence of various energy components 
(in units of hu>±) during free expansion for A = 1 (upper 
panel) and 3 (lower panel). Other parameters are Po = and 
P 2 = 1. 



energy, it always approaches zero asymptotically during 
the expansion. In the early stages, however, its behavior 
depends on trap geometry due to the anisotropic nature 
of the interaction. For small A, the initial dipole-dipole 
interaction energy is negative; during the expansion, it al- 
ways gains energy from the kinetic energy and approaches 
zero monotonically. For large A when the initial dipole- 
dipole interaction energy is positive, the condensate first 
releases its interaction energy until Pdd becomes nega- 
tive, then gains energy and finally approaches zero. 

The effects of dipolar interaction on the kinetic ener- 
gies are rather simple. We find that, independent of A, 

(r) 

the dipolar interaction always decreases P^ in while in- 
creases P^fn with time, i.e., the dipole-dipole interaction 
causes the transfer of radial kinetic energy into the ax- 
ial direction. This phenomenon is observed even when 
Po 7^ 0. This result contradicts intuition, because along 
the z-axis, the dipole-dipole interaction is attractive, one 

(z) 

would expect that the P^ would decrease with time be- 
cause the atoms are being slowed down due to the dipolar 
attraction. To resolve this puzzle, we note that the ki- 
netic energy along radial or axial directions each has two 
parts: one from the gradient of the wave function (l/2q 2 ) 
which decreases with time, and the other from the ex- 
panding gas (q 2 /2) which increases with time. This ob- 
served phenomena shows that in the radial direction the 
increase of the kinetic energy due to q 2 /2 cannot com- 
pensate for the decrease due to l/2q 2 . Similarly, we can 
understand the net increase of kinetic energy along the 
axial direction. 

For the remainder of this paper, we will focus our 
studies on 85 Rb condensate as in the JILA experiment 
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FIG. 3: The atom number dependence of the total release 
energy for a 85 Rb condensate with A = 3 (upper panel) and 
A = 6.8/17.35 (lower panel), for the s-wave scattering length 
a sc = (solid line), 0.1 (dashed line), and 0.5 (dash-dotted 
line) (as). 



JTT| , Q, where the tuning of the scattering length has 
been demonstrated in a remarkable fashion utilizing the 
Feshbach resonance. In this case, the magnetic dipole 
moment of the trapped state \F = 2, Mp = 2) is 
fi = 2/xs/3 (/ib is the Bohr magneton). We adopt the 
radial frequency loj_ = 2tt x 17.35 (Hz) as from the exper- 
iment ]IT| , |l2| and assume that the asymmetric parame- 
ter A can be adjusted. The resulting dipolar interaction 
strength is P 2 = 5.0 x KT 6 7V. 

The release energy E Ie i — Ek in + E int is the total en- 
ergy of the condensate after switching off the trapping 
potential. It's values can be strongly affected by the 
atom-atom interaction. For non-interacting atoms, the 
release energy per atom E® c} = (1 + A/2)/2, is indepen- 
dent of the atom number. For interacting atoms, it de- 
pends on the total number of atoms since both Pq and P2 
are proportional to N. If the interaction is repulsive, the 
release energy per atom always increases with N. For at- 
tractive interaction, the release energy always decreases 
(increases) with N if N is above (below) some critical 
value. Depending on the geometry of the trapping po- 
tential, the overall dipole-dipole interaction can be either 
repulsive or attractive, and the release energy also shows 
different behaviors for different trapping potentials. 

In the upper panel of Fig. |^, we present the atom 
number dependence of the release energy for A = 3. Since 
the dipole-dipole interaction is predominately repulsive 
in this case, the release energy increases with N. In the 
lower panel for A = 6.8/17.35, when P = (solid line), 
the release energy decreases with N. We find that this 
behavior holds even for a very small positive scattering 
length as shown in Fig. Q As one might have expected, 
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FIG. 4: The relative change of condensate aspect ratio due 
to the dipole interaction, (a) For the trapped ground state 
with A = 0.1, 0.5, 1, 2, 3, 4, 5, and 6. (b) The asymptotic 
condensate aspect ratio for N = 2 x 10 5 (solid line), 10 5 
(dashed line), and 4x 10 4 (dash-dotted line). Triangle markers 
show the results for N — 2 x 10° from numerically calculated 
expansions. 



only for very small values of the s-wave scattering length, 
are the effects of dipole-dipole interaction important. 

The dipolar interaction can also change the condensate 
aspect ratio. To quantify this effect, we use the relative 
difference of the (ground state) condensate aspect ratio 



A K 
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In Fig. H (a), we plot the numerical results of A K (Po = 
0, P 2 ) for various values of A. We see that the relative 
differences can become as high as 30% for a wide range 
of trap aspect ratio A. An experimental measurement 
of these differences will represent a direct detection of 
dipolar effects, although the following technical difficulty 
remains; The condensates produced in current 85 Rb ex- 
periments p"H contain rather small number of atoms that 
a in-situ direct optical imaging is very challenging. The 
free expansion as discussed in this paper, leads to larger 
condensate sizes as shown in Fig. ^, thus allowing for 
easier imaging of condensate aspect ratios. We find that 
the condensate aspect ratio changes during the expan- 
sion and eventually approaches a constant value, which 
we call the asymptotical aspect ratio. Not surprisingly, 
this asymptotical aspect ratio also depends on the dipole- 
dipole interaction. In Fig. [I| (b), we present the A de- 
pendence of the relative change of the asymptotical con- 
densate ratios for various values of N. We see that for 
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FIG. 5: The effects of dipole-dipole interaction on the stabil- 
ity coefficient k for a 85 Rb condensate with scattering length 
a sc = — 0.5ao (dashed line), — l.Oao (dash-dotted line), and 
— 3.0ao (dotted line). The solid line indicates the stability 
coefficient when dipole-dipole interaction is not included. 



carefully chosen parameters A k (Pq = 0,P-2) can become 
large enough to be observed experimentally. 

Finally, as already studied extensively before, the par- 



tially attractive nature of the dipole-dipole interaction 
can destabilize the condensate ground state in traps 
with small values of A (< 5.2) f|, g|. The stability 
coefficient k = iV cr |a sc |/aho (with aho = \Jh/muj and 
us = {oj x oj y uj z ) 1 / 3 ') is frequently used to measure the sta- 
bility of a condensate JTjJ. In Fig. [|, we plot the A de- 
pendence of k for various scattering lengths as obtained 
from variational study. We see that for a small negative 
scattering length, the dipole-dipole interaction changes 
the k value dramatically, another signature for directly 
detecting the dipolar interaction. 

In conclusion, we have studied the free expansion of 
a dipolar condensate. We have shown that the weak 
dipole-dipole interaction may become detectable through 
several observations, of which the measurement of the 
condensate aspect ratio after an expansion period seems 
the most promising. We also discussed briefly the effects 
of the dipolar interaction on destabilizing the condensate 
when the nominal s-wave scattering is tuned close to zero, 
a situation close to the recent experimental measurement 
of k in 85 Rb condensates Q. 
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